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Abstract 

In a magnetic substance the gap in the Raman spectrum, A/j, and the neu- 
tron scattering gap, Ag, are related by Ar/A^ w 2 if the the magnetic 
excitations (magnons) are only weakly interacting. But for CuGeOs the 
experimentally observed ratio is of the order A/j/A^ ~ 1.49 — 1.78, indi- 
cating attractive magnon-magnon interactions in the quasi-lD Spin-Peierls 
compound CuGeOs. We present numerical estimates of A/f/Ag from exact 
diagonalization studies for finite chains and find agreement with experiment 
for intermediate values of the frustration parameter a. An analysis of the 
numerical Raman intensity leads us to postulate a continuum of two-magnon 
bound states in the Spin-Peierls phase. We discuss in detail the numerical 
method used, the dependence of the results on the model parameters and a 
novel matrix-element effect due to the dimerization of the Raman-operator in 
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the Spin-Peierls phase. 
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The Spin-Peierls compound CuGeOa has lately been studied intensively and found 
to exhibit well defined magnetic excitations (magnons) in the dimerized phase 0. These 
magnons are found to be seperated by a gap from the continuum of two-spinon excita- 
tions predicted for the Heisenberg- chain IQ and were recently observed in KCuFs [^] and 
in CuGeOs ||^. In this context it was realized [0-0] that the magnons in the dimerized 
Heisenberg chain can be regarded as two-spinon bound-states. While spinous are essentially 
free in a homogeneous spin chain they interact strongly in a (gapped) dimerized spin chain 
where magnons are well defined excitations with dispersion uoq contributing a delta-function 
~ 5{ijj — ujq) to the dynamical structure factor, S{q,uj). 

It is then all but natural to investigate the interactions of two magnons in a dimerized 
spin-chain. Here we will present numerical and experimental evidence that magnons do 
strongly interact in dimerized spin chains leading to a continuum of two-magnon bound 
states. For this purpuse we will present data from exact diagonalization of chains with up 
to Ng = 28 sites and experimental Raman spectra for CuGeOs. We will, in particular, 
investigate the gap observed in S{q,uj) and the gap observed in the two-magnon 
Raman spectrum Ir{uj). We find generally that A^/A/j < 2, indicating strong magnon- 
magnon interactions in ID dimerized spin systems. 

As the minimal model for magnetic excitations in CuGeOs one can consider the frustrated 
ID spin Hamiltonian 

H = jY,[ii + si-iy) s, ■ s,,+i + « s,, ■ s,+2] , (1) 



where S is the dimerization parameter that vanishes above Tsp [10,11 1. The special geometry 



T^ , p!3[[ of the superexchange path in CuGeOs along the c-axis leads to a small value of the 
exchange integral J ~ 150K and a substantial n.n.n. frustration term ~ a which competes 
with the n.n. antiferromagnetic exchange. The correct value of a suitable for CuGeOs is 
still under discussion. While Castilla et al. proposed a ~ 0.24 [|10|, a much larger value 
a ~ 0.35 was proposed by Riera and Dobry [|TI|] and, recently, by Brenig et al. [Q. The 



interchain couplings have been estimated to be small, Jh ~ O.IJ and Ja ~ — O.OIJ for the 



interchain exchange constants along a- and b- directions, respectively 0. 

The phase diagram of H in Eq. ([1|) has been calculated using the density-matrix 
renormalizat ion-group method [0. For 5 = and a < Oc ~ 0.2411, the ground state is 
gapless and renormalizes to the Heisenberg fixed point. For a = 0.5 and 6 = 0, the ground 
state is given by a valence-bond state with a gap of order J/4 induced by frustration. 

An experimental method particularly suited for the study of magnetic excitations in an 
antiferromagnet is two-magnon Raman scattering. For CuGeOs, the Raman operator in Aig 
symmetry [jl6l is proportional Jl^ to 



i7^ = ^(l + 7(-ir)S.-S,+i . (2) 

i 

In the homogeneous state (5 = 7 = 0) the Raman operator commutes with the Heisenberg 
Hamiltonian for the case a = and there would be no Raman scattering p!7| , p!8| . However 
when a 7^ 0, the model (|ip leads to magnetic Raman scattering ~ due to the presence 
of competing interactions which can be observed experimentally. Note the presence of the 
factor 7 in Eq. @ which appears for T < Tsp because the exchange integral is sensitive to 
the inter-ionic distance. 

We have exactly diagonalized Eq. (|l]) for chains with up to 28 sites by a generalized 
Lanczos method and evaluated the Raman spectral weight at zero temperature, 

In{u) = -hm{0\Hn ^. ]^ ^ M O), (3) 
TT 00 + te — [H — ho) 

where Eq is the ground state energy, H the Hamiltonian given by Eq. (|1]) and e 0+. We 

have also calculated the dynamical structure factor 

S{q,u;) = J2\{n\S^JO)f6{uj-{E^-Eo)) (4) 

n 

where 5| = N-^^^ e£ 

I exp[iql]S[ and |?7,), En are the Eigenstates and Eigenenergies of the 
spin chain, respectively. 

We have evaluated S{q, u) for chains with up to Ns = 24 sites using an approximate 
scheme for the determination of the low lying excitation energies En — Eq and the corre- 
sponding transition probabilities Wn{q) = K^|'S'g|0)p [|19|. Using a recursion algorithm a set 
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of orthogonal states is built starting with 5*^10). Coefficients occuring in this procedure form 
a tridiagonal matrix whose eigenvalues and eigenstates determine the excitation energies and 
transition probabilities P,p!9|. 

For a numerical evaluation of Eq. we have used the kernel polynomial approximation. 
Since the advantages of this method have been realized only recently we give here a brief 
account. We start by rescaling the Hamiltonian by if = cX + d such that the eigenvalues of 
the rescaled Hamiltonian X are in the interval [—1, 1]. Similary we define a rescaled energy 
and frequency by Eq = cxq + d and ui = cx + d and expand Ir{x) in terms of Tschebycheff 
polynomials, Ti{x): 



= /^ ^9. + ^o), (5) 



Vl- x^ 



where the number of polynomials retained, Np, determines the accuracy of the approximation 
which becomes exact in the limit Np oo. The expansion coefficients a; are determined 
using the orthogonality relations for Tschebeycheff polynomials to be 

TT 

which can be evaluated recursively via the formula T;+i(x) = 2xTi{x) — T/_i(x). The advan- 
tage of an expansion in orthogonal polynomials is its numerical stability. It is indeed possible 
to evaluate several thousands of ai recursively without encountering numerical instabilities. 
Often we will use only a limited number Np = 100 for comparison with experimental data. 

A truncated expansion in orthogonal polynomials will, in general, lead to unwelcome 
Gibbs oscillations for any finite Np < oo. These Gibbs oscillations have been studied care- 



fully in the past |20| and can be suppressed efficiently and reliable by the replacement 



O'l — * ciig{zi) with g{zi) = [sin^TTZi)/ {irzi)]'^ and zi = l/{Np + 1) [^]. Note that this replace- 
ment, which we use throughout this paper, still satisfies the correct limit Np ^ oo. 

A finite A'p < cxd in does broaden the delta-poles in a finite system calculation. The 
dependence of Ir{uj) on Np is illustrated in Fig. |I| for a chain with Ng = 28 sites, a = 0.24 
and 5 = = 7. Ill the inset of Fig. |l| we compare the broading of a single pole with the kernel 



polynomial approximation {Np = 100) with a Lorentzian e/{{uj — cjj)^ + e^) and e = 0.025. 
Note the absence of the high-energy tails in the kernel polynomial approximation. 

We have measured the Raman intensity on single crystal CuGeOa using the excitation 
line A = 514.5-nm of an Ar-laser with a laser power of 2.7mW. We ensured that the incident 
radiation does not increase the temperature of the sample by more than 1.5 K. We used a 
DILOR-XY spectrometer and a nitrogen cooled CCD (back illuminated) as a detector in a 
quasi backscattering geometry with the polarization of incident and scattered light parallel 
to the c-axis and the Cu-0 chains, respectively. 

In Fig. 1^, we present the data for the two-magnon Raman continuum in the homogeneous 
state at T = 20 K. Phonon lines [p2| , p3| at 184cm~^ and at 330cm~^ are subtracted from 
the experimental data (squares). The experimental Raman spectrum presented in Fig. Q is 



in agreement with other Raman studies on CuGeOa p2[. We have included in Fig. y the 
numerical results for Ir{uj) obtained for chains with A^^ = 24 (dashed lines) and Ns = 28 
(solid lines) sites for 5 = = 7 ^-^id Np = 100. Note that the finite-size effects are quite 
small. We show in Fig. |^ data for two parameter sets, namely a = 0.24, J = 150K and 
a = 0.35, J = 159K. We note that a = 0.35, which is favoured by fits to the susceptibility 



| rT| and to the specific heat does not agree well with the Raman spectrum. Similar 



results have been obtained previously with a solitonic mean-field approach to the frustrated 
Heisenberg chain [[1^ , ^^] . 

In the dimerized phase 5 7^ we have found that the numerically obtained Raman 
spectrum depends very much on the dimerization parameter 7 in the Raman operator as 
we illustrate in Fig. ^ for a chain with A^^ = 24 sites, 6 = 0.03, a = 0.24 and Np = 100. 
Between 7 = and 7 = 0.15 the spectrum changes qualitatively and a low-energy peak can 
be resolved for 7 = 0.12, 0.15, but not for 7 = 0. 

In order to understand the dramatic dependence of the Ir{uj) on 7 observed in Fig. 
H we have analyzed the dependence of the seven lowest poles contributing to Ir{uj) on 7, 
illustrated in the inset of Fig. | for A", = 24, a = 0.24, 6 = 0, Np = 1000 and 7 = 0. 

We start by rewriting Eq. (^ in the form 



Ir{uj) = Y.\{0\Hn\n)\'5iuj-iE^-Eo)) (6) 

n 

and noting that we can decompose the Raman operator into two parts, Hji = H' + jH". 
The weight of an excited state then becomes 



\{0\H' + ^H"\n)\'' = \m' + -fm"\\ (7) 



and for any n there is a 70 = —m'/m" for which w„ vanishes. We have analyzed the energy- 
dependence of 7o = 7o(-E'„) and found that for the five dominant poles, n = 1, 2, 4, 6, 7 
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7o ^ 5 + En, (8) 

in inverse wavelength [cm^^] for En- This dependence of w„ on En is shown in Fig. § where 
we have plotted the weights as a function of 7 — 70 (-En)- The data presented in Fig. | clearly 
indicates that the dominant contributions to Ir{uj) follow the scaling relation 



l-y=0 (1 - loiEi) 
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where is a normalization constant, which we have approximated by the constraint 
Jq" dujp{-y,Lj) = 1, with Uc = 6 J. The scaling relation Eq. (|^) constitutes, on the other 
hand, also the rescaling of the Raman intensity Ir^uj) = Ir{j,u!) with 7, such that we can 
write 

/fi(7,^)^p(7,a;)/fi(0,^). (10) 

We see that the effect of 7 is to "burn a spectral hole" in Ir{0, uj) at a characteristic frequency 
7o(ct;), in agreement with the numerical results presented in Fig. |^. 

We can actually find a complete analytic formula for /i?(7, uj) by noting that for a = 0.24, 
6 = 0.03 we can approximate Ir{0,u!) by the expression (see Fig. || and |I7|] ) 



Jh(7 = 0, ^) = Ae{u - Ar) (1 - tanh[2(^ - cuo)]) , (11) 

with the values of the parameters Ar ^ 30cm^^, ujq ~ 312cm~^ being determined by a fit 
to the numerical data (A is a normalization constant). The heavy-side function 6{uj — Ar) 
in Eq. (|Tl|) reflects the absence of two-magnon excitations below Ar in the dimerized state. 



In Fig. ^ we present the experimental Raman spectrum in the spin-Peierls phase at 
T = 5K (with phonon hnes substracted, filled squares) in comparision with the analytic 
result for 7 = 0.12(Eq. (|T0|) together with Eq. (|11]), solid line). Experimentally a two- 
magnon peak is observed at 30cm~^ and a broad continuum at higher frequencies. These 
two features are well reproduced by the analyitic result. There is, on the other hand, a new, 
probably magnetic line at 225cm^^, which is not reproduced by our theoretical study of the 
one-dimensional spin model (|l|), compare also Fig. |. This fact has led us to speculate fl^ 
that the interchain coupling Jb might become relevant for T < T^p in CuGeOa. 

Which is the physical significance of the peak at Ap ~ 30cm~^ in Fig. ^ observed both 
in experiment and in our analytical result? Neutron-scattering experiments [0,^ indicate 
a one-magnon gap A5 ~ (2.1 — 2.5)meV, i.e. ~ (16.9 — 20.2)cm~^. The experimental 
ratio 



As 



~ 1.49 - 1.78 (12) 

Exp. 



is smaller than two, the value expected for non-interacting magnons. This indicates an 
attractive interaction between magnons in ID dimerized spin chains, similar to the attraction 
between spinous. Uhrig and Schulz |^ have indeed postulated an isolated singlet bound state 
with energy Eg and 2As' > > A^ which might be a two-magnon bound state. 

A singlet bound-state with zero total momentum would be observable by Raman scat- 
tering and would show up as a delta-function contribution to the Raman intensity (H): 

Ir{uj) = A5{uj-E,) + 4(^), 

with I'{i{uj) beeing the continuum part of the Raman intensity. A finite weight A implies 
also a finite value of the relative weight of the first pole, 

Ar = lim (13) 

with Wn = I (0|iJij|n) p. In the inset of Fig. ^ we show the relative weight as a function of 
1/A^^ for chains with up to A^^ = 28 sites and a = 0.24. We find no indication for an isolated 
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bound state, i.e. for a finite value of Aj. both for S = (filled circles), as expected, and for 

5 = 0.03 (filled triangles). The data for 6 = 0.03 shown in the inset of Fig. ^ have been 
calculated for 7 = 0.12. The data for 7 = and 6 = 0.03 are very similar. 

We have therefore no indication from numerics for an isolated singlet bound state, in 
agreement with our interpretation of the 30cm~^ peak as a continuum contribution truncated 
by matrix-element effect, see Eq. (p!oD. Our numerical results show, on the other hand, 
unambigously that the lower edge of the Raman spectrum is pulled below the non-interacting 
two-magnon density of states. This fact, which is discussed in detail further below, indicates 
a continuum of two-magnon bound states with zero total momentum. This scenario is 
conceivable in view of the fact, that the two-magnon spectrum is made up out of four 
spinous 

The experimentally observed values of single- and two-magnon gap energies can be used, 
to a certain extent, to determine the value of the parameters J, 6 and a entering Eq. (p. 
The ratio Aji/ As is independent of the coupling-constant J. We have plotted Aji/ As in Fig. 
^ for various values of a as a function of 6. The data are obtained for Ng = 24, the finite-size 
corrections are very small (compare also [^). The experimental possible values 1.49-1.78 are 
indicated as the shaded region in Fig. ^. Again, we see that a larger value of a ~ 0.35 does 
not do a good job on An/ As- In using Fig. |]for determing possible combinations of a and 

6 one has to keep in mind that the interchain coupling J^, not included in our caculations, 
has an effect of order 10% on the values quoted in Fig. ^ It is interesting to note, that 
recent experiments on CuGeOs under pressure ||2^ show that the frustration parameter a 
increases with presure, as does the experimental ratio Aji/As, in accordance with the data 
presented in Fig. ||. 

In conlusion we have discussed in detail the effect of the dimerization of the Raman 
operator on the Raman spectrum in the spin-Peierls state of CuGeOs. We have not found 
any numerical evidence for a singlet bound-state and conclude that the 30cm~^ Raman line 
observed at T = 5K in CuGeOa is a two-magnon line which is carved out of a continuum 
by matrix-element effects due to the dimerization of the Raman operator, 7. We have then 

9 



compared the observed ratio of the Raman gap to the neutron-scattering gap and compared 
with numerical results. We conclude that this ratio indicates an intermediate value for the 
frustration parameter a. 
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FIGURES 

FIG. 1. Dependence of the spectrum obtained by the kernel polyonmial approximation on the 
number of Tschebeycheff polynomials, Np = 100, 300, 800, retained. The data is for a chain with 
Ng = 28 sites and 5 = = 7, a = 24. For Np — > 00 the poles evolve into delta- funtions. Insert: 
Comparison between a pole located at x = 0.5 broaded by a Lorentzian e/{{x — 0.5)^ + e^) (dashed 
line), with e = 0.025 and an orthogonal polynomial approximation with Np = 100 Tschebycheff 
polynomials. 

FIG. 2. The magnetic contribution to the experimental Raman spectrum for CuGeOs at 
T = 20K (filled squares) and data obtained for chains with Ns = 24,28 sites (dashed/solid line) 
for a = 0.24 and a = 0.35 respectively and Np = 100 and 5 = = 7. 

FIG. 3. Dependence of the spectrum in the dimerized phase, S = 0.03, a = 0.24, 
Np = 100 for a chain with Ns = 24 sites on the dimerization of the Raman-operator, 
7 = 0, 0.03, 0.06, 0.09, 0.12, 0.15. Note the large matrix-element effects. 

FIG. 4. The rescaled weight w{En) of the lowest seven poles in the Raman spectrum with 
pole Energy, En (n = 1, . . . , 7) as a function of the rescaled dimerization of the Raman operator, 
7 — jo{En), see Eq. (^). The data is for a chain with Ng = 24 sites and 6 = 0.03, a = 0.24. The 
weight is rescaled by 7o(£'i) = 6 + KEi, with n ~ 4/3000cm. We observed that the rescaled weights 
of the dominant poles 1, 2, 4, 6, 7 fall all on a universal curve (the data for pole number 4 has been 
omitted in order to avoid overcrowding). 

FIG. 5. Comparision between the experimental Raman spectrum in the dimerized phase of 
CuGeOs at T = 5K (filled squares) and the analytic expression ((p!o|), solid line) Insert: The 
relative weight (Eq. (|l3|)) of the lowest pole contribution to Ir{uj) as a function of l/A^f for 
Ns = 12, 16,20,24,28. The data is for a = 0.24, <5 = (filled circles) and 5 = 0.03 and 7 = 0.12 
(filled triangles). 
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FIG. 6. The calulated ratios of the Ar/As from chains with Ng = 24 sites for various a and 
5. The shaded region indicates the experimentally observed values for CuGeOa. 
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